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Enzymatic  reactions  in  biochemistry  are  far  more  complex  than  the  celebrated  MM  scheme,  but  the  observed  turnover  rate 
often  obeys  the  hyperbolic  dependence  on  the  substrate  concentration,  a  relation  established  almost  a  century  ago  for  the 
simple  MM  mechanism.  Recent  single  molecule  studies  of  enzymatic  reactions  reveal  that  proteins  fluctuate  over  a  wide  range 
of  time  scales,  but  turnover  reactions  of  fluctuating  enzymes  are  often  found  to  obey  the  MM  relation,  which  was  originally 
derived  for  enzymes  with  a  unique  conformation.  To  resolve  the  longstanding  problem,  we  apply  the  flux  balance  method  to 
predict  the  functional  form  of  the  substrate  dependence  in  the  mean  turnover  time  of  complex  enzymatic  reactions  and  identify 
detailed  balance  (i.e.  the  lack  of  unbalanced  conformational  current)  as  a  sufficient  condition  for  the  MM  equation  to  describe 
the  substrate  concentration  dependence  of  the  turnover  rate  in  an  enzymatic  network,  thus  establishing  a  connection  to  a 
fundamental  property  of  non-equilibrium  thermodynamics. 

When  the  DB  condition  for  the  conformational  states  is  violated,  the  MM  functional  form  will  break  down,  and  the  turnover  rate 
will  depend  on  the  unbalanced  current.  Though  explicit  rate  expressions  are  difficult  to  obtain,  we  established  a  functional  form 
for  the  mean  turnover  time  as  MM  terms  plus  non-MM  correction  terms  arising  from  the  presence  of  the  non-equilibrium  current 
in  the  cyclic  loop.  In  fact,  we  can  rigorously  prove  the  one-to-one  correspondence  between  the  unbalanced  currents  and  the 
independent  hyperbolic  terms  and  use  this  result  to  demonstrate  a  broad  range  of  cooperative  behaviors  in  non-equilibrium 
enzymes.  Using  this  generalized  MM  rate  expression,  we  were  able  to  establish  a  phase  diagram  of  different  kinetic  behaviors 
in  enzyme  networks,  including  positive  cooperativity,  substrate  inhibition,  bi-stability,  and  oscillations.  This  finding  helps  analyze 
recent  single  molecule  studies  of  enzymatic  networks  and  can  be  applied  to  other  external  variables,  such  as  force-dependence 
in  molecular  motors  and  voltage-dependence  in  ion  channels. 

A  major  part  of  work  has  been  published  in  a  JPC  special  issue  and  in  “Advances  in  Chemical  Physics  on  Single  Molecule 
Biophysics:  Experiments  and  Theories”,  a  volume  in  the  Advances  in  Chemical  Physics  series  published  by  John  Wiley 
publishers.  Two  more  manuscripts  are  now  ready  for  submission:  one  on  the  study  of  detailed  balance  violation,  and  the  other 
on  the  general  functional  form  of  substrate  dependence  in  reversible  enzymatic  reactions. 


b)  Counting  statistics  of  single  molecule/cell  measurements 

Fluctuations  in  reactivity  of  biological  molecules  are  ubiquitous,  and  their  effects  on  single  molecule  reaction  trajectories 
become  observable  with  the  advent  of  modern  single  molecule  experimental  techniques.  The  emergent  challenge  in  this  field  is 
to  extract  robust  quantitative  information  from  the  noisy  time  series  of  single  molecule  reaction  systems.  In  this  project,  we  focus 
on  a  universal  kinetic  transition  in  reaction  event  counting  statistics.  By  analyzing  the  Poisson  indicator  in  the  reaction  event 
counting  statistics,  one  can  obtain  quantitative  information  about  the  magnitude  and  the  relaxation  dynamics  of  reactivity 
fluctuations  of  biopolymer  molecules  without  presuming  a  particular  model  for  the  conformational  dynamics  and  reaction  rate  of 
a  biopolymer. 

Related  to  single  molecule  measurements,  the  single  cell  essay  technique  poses  similar  challenges  in  data  analysis  and  defines 
a  new  area  of  statistical  study.  Specifically,  the  single  cell  microarray  measurement  collect  large  amount  of  data  on  protein 
expression,  which  is  regulated  by  protein  kinases  of  Plasmodium  falciparum.  We  have  developed  a  theoretical  formulation  to 
de-convolute  and  synchronize  microarray  data  for  better  characterization  of  underlying  genetic  control  involved  in  the  malaria 
infection  cycle. 

Since  the  last  report  in  201 1 ,  this  project  has  become  a  major  direction  in  my  group.  Two  more  papers  have  since  been 
published  on  single  cell  counting  statistics.  In  2013,  we  published  a  JPC  paper  jointly  with  a  Los  Alamos  group  on  ‘universality 
of  Poisson  indicator  and  Fano  factor  of  transport  event  statistics’.  This  paper  will  be  followed  by  another  manuscript  on  the 
generic  scheme  for  enzymatic  reactions.  In  these  two  papers,  we  systematically  examine  the  functional  form  of  substrate 
dependence. 


c)  Joint  theoretical-experimental  study  of  cytoadhesion  and  protein-receptor  binding  kinetics 


The  project  is  a  part  of  the  recent  efforts  in  my  group  to  develop  a  joint  theoretical-experimental  program  in  biophysics,  including 
cell  mechanics  and  adhesion.  Cell  adhesion  plays  a  very  important  role  in  multi-cellular  organisms,  helping  to  maintain  their 
integrity  and  homeostasis.  This  complex  process  involves  many  different  types  of  adhesion  proteins,  and  the  synergetic 
behavior  of  these  proteins  during  cell  adhesion  has  been  frequently  observed  in  experiments.  A  well-known  example  is  the 
cooperation  of  rolling  and  stationary  adhesion  proteins  during  the  leukocytes  extravasation.  Despite  the  fact  that  such 
cooperation  is  vital  for  proper  functioning  of  the  immune  response,  its  origin  is  still  insufficiently  understood.  In  this  study,  we 
constructed  a  simple  analytic  model  for  the  interaction  between  leukocytes  and  blood  vessel  walls  in  shear  flow  and  predict  the 
bistability  of  cell  adhesion.  This  bistabilty  results  from  the  competition  between  two  kinetic  processes  in  the  cell-wall  contact 
area  -  bonds  formation  and  rupture.  The  theoretical  model  shows  how  physicochemical  properties  of  adhesion  proteins 
influence  the  adhesion  bistability.  Based  on  these  results,  we  propose  an  explanation  of  the  synergy  between  rolling  and 
stationary  adhesion  proteins  for  effective  cell  adherence  to  the  blood  vessel  wall  in  shear  flow  and  provide  possible 
interpretations  of  several  cytoadhesion  experiments. 

The  theoretical  model  on  cell-adhesion  was  published  in  Biophysical  Journal,  and  a  joint  theory-experiment  manuscript  was 
published  in  PLOS  One.  Currently,  our  effort  is  directed  to  the  detailed  study  of  the  binding  kinetics  of  cell-substrates,  relating 
cytoadhesion  to  protein-receptor  interaction.  This  effort  leads  to  another  manuscript  under  consideration  in  Biophysical  Journal. 


d)  Coherent  energy  transfer  in  photosynthetic  light-harvesting  systems 

In  photosynthesis,  excitation  energy  is  transferred  from  pigments  of  the  antenna  to  the  reaction  centers  with  remarkable 
efficiency.  Recent  experiments  reveal  that  photosynthetic  systems  and  their  surrounding  protein  environments  interact 
coherently  to  optimize  transfer  pathways  and  maximize  efficiency.  The  high  efficiency  in  natural  systems  has  inspired  scientists 
to  design  artificial  systems  for  efficient  conversion  of  solar  energy  into  useful  energy  forms.  These  experimental  results  have 
motivated  two  theoretical  questions:  how  to  optimize  the  efficiency  of  energy  transfer  systems  in  dissipative  environments  and 
how  to  understand  the  optimization  conditions  in  terms  of  network  kinetics.  In  several  recent  manuscripts,  we  developed  two 
theoretical  approaches  to  address  these  questions:  generic  mechanism  of  optimal  energy  transfer  and  systematic  expansion  for 
mapping  exciton  dynamics  to  network  kinetics. 

The  optimization  of  excitation  energy  transfer  depends  on  initial  preparation,  static  disorder,  and  the  type  of  dynamic  noise.  To 
unify  these  different  physical  conditions  into  a  single  conceptual  framework,  we  introduced  the  concept  of  orthogonal  subspace, 
which  is  the  part  of  exciton  space  that  is  decoupled  from  the  trap  state.  In  the  coherent  regime,  the  energy  transfer  from  an 
orthogonal  exciton  subspace  is  prohibited  and  needs  assistance  of  dynamic  and  static  disorders.  In  the  incoherent  regime, 
excitation  energy  can  be  localized  due  to  the  dissipative  effect  from  the  surrounding  environment.  The  competition  between 
these  two  mechanisms  thus  leads  to  an  optimization  at  an  intermediate  level  of  noise.  The  orthogonal  subspace  is  a  general 
concept  for  understanding  optimal  efficiency  and  can  be  found  in  systems  with  topological  symmetry,  limited  quantum  coherent 
strength,  and  large  system  size. 

To  interpret  the  analytical  solutions  for  multi-level  systems  and  the  numerical  solutions  for  realistic  photosynthetic  systems,  we 
constructed  a  systematic  expansion  to  map  exciton  dynamics  to  more  intuitive  network  kinetics.  The  expansion  defines  an 
effective  hopping  rate  as  the  leading  order  picture  and  non-local  kinetic  couplings  as  the  quantum  correction,  leading  to  a 
rigorous  separation  of  thermal  hopping  and  coherent  tunneling,  and  is  a  useful  tool  for  visualizing  pathway  connectivity  and 
interference  in  quantum  networks.  The  mapping  procedure  allows  us  to  apply  the  analytical  tools  developed  for  chemical 
networks  to  study  photosynthetic  energy  transfer.  We  are  currently  working  on  a  general  kinetic  expansion  of  quantum 
dissipation,  which  turns  out  to  be  equivalent  to  a  generalization  of  NIBA  (non-interacting  blip  approximation)  for  the  spin-boson 
model  to  NICA  (non-interacting  cluster  approximation)  for  multi-level  exciton  systems. 

Late  Professor  Silbey  was  a  leader  in  this  research  direction.  Together,  we  published  over  10  papers  on  this  subject,  including 
papers  on  PRL,  Nature  Communications,  PNAS,  and  JPC  letters.  Also,  we  contributed  to  chapters  in  the  recent  book:  ‘Quantum 
Effects  in  Biological  Systems'.  This  project  was  sponsored  jointly  by  this  ARO  grant  (joint),  NSF  grant  (individual),  and  a 
DARPA  grant  (joint). 
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